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Abstract. We present molecular line and submillimeter dust continuum observations of the Lynds870 cloud in 
the vicinity of IRAS 20231+3440. Two submillimeter cores, SMM1 and SMM2, are identified mapping the 870 /xm 
dust continuum and ammonia (NH3) emission. The total molecular mass is ~70 - IIOMq. The northern core is 
warmer and denser than the southern one. Molecular outflows are discovered in both cores. In the northern one a 
significant amount of low velocity (1.3-2.8 km s _1 ) outflowing gas is found, that is hidden in the relatively broad 
CO lines but that is revealed by the narrower HCO + spectra. While IRAS 20231+3440 is most likely the exciting 
star of the northern outflow, the driving source of the southern outflow is not detected by infrared surveys and must 
be deeply embedded in the cloud core. Large scale (~0.2pc) infall motion is indicated by blue asymmetric profiles 
observed in the HCO + J = 3-2 spectra. Red K s band YSO candidates revealed by the 2MASS survey indicate 
ongoing star formation throughout the cloud. The calculated masses and the measured degree of turbulence are 
also reminiscent of clouds forming groups of stars. The excitation of the molecular lines, molecular abundances, 
and outflow properties are discussed. It is concluded that IRAS 20231+3440 is a Class I object, while the southern 
core most likely contains a Class source. 
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1. Introduction 

Stars form from collapsing dense cores of molecular clouds. 
Protostars arise as a result of the collapse and evolve (if 
their mass does not exceed a few solar masses) from Class 
to Class I, II and III objects (Andre et al. |2000| ). The 
earliest stages of star formation are observationally char- 
acterized by an association with dense molecular cores, 
masers and molecular outflows. The dust and the gas of 
dense cores can be traced by submillimeter continuum 
and spectral line emission from tracers like NH3, CS, and 
HCO + . Molecular outflows, usually traced in lines of the 
CO molecule with broad-line wing emission, have received 
considerable attention over the past two decades and over 
250 outflow sources have been cataloged so far (Fukui et 
al. 1993; Wu et al. 1996). Our current understanding of 



dense core is, however, heavily influenced by a few well- 
studied examples, and may be far from accurate. 

Recently an IRAS-selected CO J = 1-0 survey has 
been carried out at the 13.7m mm-wave telescope of the 
Purple Mountain Observatory at Delingha (Yang et al. 
1999| , |2002j ; Jiang et al. |2000p , A number of sources with 



outflow properties and its interaction with the parental 
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Telescope (HHT) and the MPIfR 100-m telescope at Effelsberg. 
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Observatory on behalf of Steward Observatory and the Max- 
Planck-Institut fur Radioastronomie. The 100-m telescope is 
operated by the Max-Planck-Institut fur Radioastronomie. 



strong CO emission and/or outflowing gas have been se- 
lected from the survey. In this work, we present detailed 
observations of one of the sources, IRAS 20231+3440 
(hereafter 120231). 

120231 (a 200 o.o=20' , 25 m 07 s .O, 5 2 ooo.u=34°50'05") is 
located within Lynds870. On Palomar Observatory Sky 
Survey (POSS) prints, there is no optical counterpart at 
the position of 1 20231. High extinction can be traced over 
a region of 15' x 15'. The infrared colour indices of 1 20231 
show that the source is deeply embedded (Emerson 1988). 
At the same time the IRAS colors satisfy the criteria in- 



troduced by Wood & Churchwcll (1989) for compact HII 
regions. There were some water maser detections reported 



(Palla et al. 1993; Brand et al. 1994), which indicate an 
early stage of young stellar object (YSO) formation, but 
methanol (van der Walt et al. 1996) and OH masers were 



not detected (Slysh et al. 1997). CS 2-1 observations give 
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fLSR^5.4kms _ (Bronfman ct al. 1996), which indicates 
a dynamical distance of 3.7 kpc (far) or 1.0 kpc (near dis- 
tance). The optical extinction seen on POSS prints sug- 
gests that the nearer distance is more reasonable and we 
will take 1.0 kpc as the distance to the source. Since all 
mapping observations were carried out w.r.t 120231, we 
refer in the following to this as the reference position (off- 
set (0", 0")). 

In this paper, we present molecular line and 870 fim 
continuum observations to trace the dense molecular 
gas, the outflows and the submillimeter dust continuum 
around 120231. In §||, we describe our observations, and 
in and §|I] we derive physical parameters from our 
dust continuum and molecular line data. In § |, results 
from NIR photometry are studied and in §p, we dis- 
cuss the global properties of star formation in the region. 
Conclusions are drawn in 

2. Observations 

2.1. NH 3 and H 2 at Effelsberg 100-m 

The NH3 and the 22 GHz H2O maser observations were 
obtained in July 1999, April and July 2000, using the 
Effelsberg 100-m telescope of the Max-Plank Institute fur 
Radioastronomie (MPIfR) equipped with a maser receiver 
in 1999 and a dual channel K-band HEMT receiver in 
2000. The typical system temperature was about 150 K - 
200 K on a main beam brightness temperature scale; the 
half-power beam width was ~40". The data were recorded 
with a 8192-channel autocorrelator. The K-band HEMT 
receiver allowed us to measure the (J, K) = (1,1) to (4,4) 
NH3 inversion lines simultaneously. The observed spectral 
resolution was 0.24 and 0.98kms _1 for NH3 observations 
in 1999 and 2000, respectively, and about 0.13 km s -1 for 
the H2O maser observations. The line intensities were cal- 
ibrated relative to the continuum of NGC 7027, for which 
we assume a flux density of 5.86 Jy, corresponding to a 
main beam brightness temperature of 8.1 K. Pointing was 
checked every hour on nearby continuum sources and was 
found to be stable within 10". The calibration accuracy is 
estimated to be ±20%. 

2.2. mm and sub-mm line observations at HHT 10m 

Molecular line emission of CO, CS, HCO + and some of 
their rare isotopes was observed with the 10-m Hcinrich 
Hertz Telescope (HHT; see Baars et al. |1999Q on Mt. 
Graham in Southern Arizona during three sessions in 
April 1999, January 2000 and May 2000. The receivers 
were sensitive to both sidebands. A dual channel 345 GHz 
SIS receiver was used for the J = 3-2 transitions of 12 CO 
and its isotopes, while all the mm-wave lines were ob- 
served with a single channel 230 GHz SIS receiver. For 
the 12 CO mapping we used the spectral line On-The-Fly 
(OTF) method, in which we scanned the telescope at a 
rate of about 1" per second. The off position was located 
at Aa = 1° west of 1 20231 and was free of CO emission. 



The most central region (~l'-2') around 120231 was also 
mapped in 13 CO, C 18 0, C 17 0, HCO+, H 13 CO + and CS 
with standard position switching (off position at Aa = 
20')- The J = 5-4 lines of three rare CS isotopes, 12 C 34 S, 
12 C 33 S and 13 C 32 S, were obtained towards the central po- 
sition of 120231. 

For the 12 CO observations we used two 2048-channel 
1 GHz Acousto-Optical Spectrometers (AOS) as backends, 
with a mean channel spacing of 917 KHz corresponding to 
-0.78 and 1.19 kms -1 at 345 and 230GHz, respectively. A 
256-channel filter bank with a channel spacing of 250 kHz, 
providing higher spectral resolution, was used simultane- 
ously for other molecular lines. The beam widths were 
about 22 and 33" for 345 and 230 GHz respectively. 

All results displayed are given on a main beam bright- 
ness temperature scale (Tmb), whi ch is r elated to via 
T MB = ZI (F e $/B eS ) (see Downes |l989| ). Main beam ef- 
ficiencies, B c g, were 0.5 and 0.78 at 345 and 230 GHz, 
as obtained from measurements of Saturn. Forward hemi- 
sphere efficiencies, F c g, were 0.9 and 0.95, respectively. 
For 12 CO, xCyg was used for both pointing and absolute 
line intensity calibration. The rms pointing uncertainty 
was found to be ~5". Line intensities were found to be 
accurate within 20% when compared to the J = 3-2 line 



of CO of Sta nck ct al. (1995) and the J = 2-1 line of Loup 
et al. ( 1993 ). For other lines, Orion-KL was observed as a 
second ary ca librator and the lin e surv ey results by Sutton 
et al. ( |1985| ) and B lake et al. ( |l986|) near 230 GHz and 
Schilke et al. ( 1997 ) near 345 GHz were used for compari- 
son. Again calibration errors should not exceed 20%. 

2.3. 870/im continuum observations 

870 /xm dust continuum emission was imaged with the 
HHT 19 channel bolometer array. The array covers a 
hexagonal 200" field of view. The beam width of the HHT 
is about 22" at 870 ^m. Maps were taken using the OTF 
method with a beam throw of 120" and a scanning speed 
of 8" per second. Five OTF maps, with three coverages 
(400" x 300" in size) centered at 120231 and another two 
coverages (320" x 260") centered at (0",-180") relative to 
120231 were finally combined by mosaicing. The average 
r.m.s noise level is about 90mJy in the combined map. All 
maps were calibrated using 870 fim skydips and the stan- 
dard calibrators Uranus and G 34.3+0.2, and a peak flux 
density of 60 Jy beam -1 was assumed for the latter one. 
At 870 fim, the typical atmospheric opacity was about 0.4 
during the observations. The calibration accuracy should 
be better than ±20%. 

3. Dust Continuum Emission 

A 870 urn. continuum image of L 870 is presented in Fig. 
The map reveals two submillimeter cores (SMM), the 
northern one, SMM1, peaked at (7", -7") with a peak flux 
density of 1.67 Jy beam -1 , directly associated with 1 20231 
and the southern one, SMM2, peaked at (12" -188") with 
a peak flux density of 0.60 Jy beam -1 , not associated with 
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Fig. 1. HHT 870 /im continuum image overlaid on a 
2MASS-X S image. Contour levels are 15% (~3cr) to 95% 
of the peak flux (~1.7 Jy) by steps of 10%. The position 
of the IRAS point source 1 20231 is indicated by the white 
cross and the ellipse is the IRAS la positional uncertainty. 
Also shown are the 12 red K s sources IRS1, 2 to 12, which 
are discussed in § |[ 



any IRAS point source. The northern core region is elon- 
gated along a position angle of P. A. ~ -20° with a decon- 
volved half-power size of 17" by 50" (0.08 pc by 0.25 pc), 
while the southern core region looks more irregular. Also 
shown in the map is a filament southwest of 1 20231. This 
is however only a ~3cr detection, and deeper images are 
needed to confirm its existence. The integrated flux densi- 
ties of the northern, the southern and the filament regions 
are 6.5±1.1, 1.7±0.5 and 0.7±0.2 Jy, respectively. 

The masses of the clouds can be estimated via 
M=F u LP/K u B u (T d ), where F v is the flux density D de- 
notes the distance towards the cloud, -B„(Td) is the Plank 
function for temperature Td, and k„ is the dust opac- 
ity per unit gas and dust mass. Various values of k v 
for different regions can be found in the literature (e.g. 
Draine & Lee 11984; Hildebrand |1983l Preibisch et al. 



1993; Ossenkopf & Henning 1994; Kriigel & Siebenmorgen 
1994). We adopt here the value K87o^im = 0.01cm 2 g _1 , 



using the Hildebrand (1983) assumptions, with {3=2 (/?: 
exponent in the frequency dependence of the dust grain 
emissivity). For the southern dust core, we apply a dust 
temperature of 12 K, a typical value for star forming cores 



(Williams et al. |1999| Bacmann et al. |200Cj ). For the north- 
ern core, we take a higher dust temperature of 30 K, as in- 
dicated by the color temperature of ~ 33 K derived from 
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Fig. 2. Effelsberg NH3 map (contours) of the integrated 
(1,1) main group of lines overlaid on the 870/im continuum 
image (grey). Contour levels are 20% to 90% of the peak 
value (9.5Kkms _1 ) with a spacing of 10%. NH 3 (1,1) to 
(4,4) spectra taken at the 1 20231 position are attached. 
Note that the spectral resolution is 0.24 km s^ 1 for the 
(1,1) and (2,2) lines and 0.98 km s -1 for the (3,3) and (4,4) 
lines. For the white cross and the error ellipse, see Fig|l]. 



60 and 100/im dust emission (Henning et al. 1990 ). Total 
masses traced by the dust are then 33.9±6.3, 40.0±11.8 
and 16.5±4.7M0 for the northern core, the southern core 
and the filament, respectively. Note that the masses de- 
rived in this way can be linearly scaled to another value 
of K870/im if desired; moreover masses are sensitive to the 
temperature and a higher temperature would decrease the 
calculated mass inversely proportional. 

The beam averaged (22") peak column density, de- 
rived from the peak flux density of 1.67 Jy beam -1 , is 
9.6xl0 22 cm -2 , which gives a corresponding peak extinc- 
tion of A 
al 



■ 100 (iV(H 2 ) = A v x 9.5 x 10 20 cm- 2 , Bohlin et 
1978). The beam averaged peak column density for a 



different beam size can be calculated by smoothing the im- 
age to the corresponding beam size, assuming a Gaussian 
shape of the beam. Column densities of 7.5, 7.0, 6.3 and 
5.5xl0 22 cm~ 2 are derived for beam sizes of 29, 32, 35 and 
40", respectively To estimate molecular abundances as- 
suming they trace the same volume as the 870 /*m dust 
continuum emission, these values will be used in the fol- 
lowing. 
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Table 1. Observed molecular line parameters* 



Line 


Frequency 
(MHz) 


HPBW 

(arcsec) 


Tmb 
(K) 


VhSR 

(kms -1 ) 


AVint 

(kms" 1 ) 


J TMBtfo 

(Kkms" 1 ) 


T 


Offset 
(",") 


12 C le O(3-2) 


345795.975 


22 












«A 8 


.Zo ) 




(0,0) 


12 C 16 0(2-1) 


230537.990 


33 












89 8 
oZ.o 


^9^ 




(0,0) 


13 C 16 0(3-2) 


330587.957 


23 


5.37 (.12) 


D.UU 




'i A K 




in Q 
ly .o 


1 9"! 


3.78 (1.54) 


(0,0) 


13 C 16 0(2-1) 


220398.686 


35 


5.76 (.12) 


f\ i n 

D. 1U 




O.Oo 




99 « 


1 CA 


3.84 (1.00) 


(0,0) 


12 C ls O(3-2) 


329330.570 


23 


2.20 (.15) 


6.17 




3.19 




8.03 


.30) 


0.50 (0.21) 


(0,0) 


12 C ls O(2-l) 


219560.360 


35 


3.31 (.14) 


6.24 


(.03) 


2.60 


(.06) 


9.26 


.20) 


0.51 (0.13) 


(0,0) 


12 C 17 0(2-1) 


224714.368 


35 


0.45 (.06) 


6.60 


(.15) 


1.65 


(.45) 


1.41 


.15) 




(0,0) 


12 C 32 S(5-4) 


244935.643 


32 


1.67 (.14) 


6.26 


( 05) 


3.23 


f 13) 


5.77 


15) 


1.81 (0.84) 


(0,0) 








0.51 (.13) 


6.08 


( 10) 


2.72 


( 241 


1.49 


11) 




(0,-180) 


12 C 34 S(5-4) 


241016.194 


32 


0.13 (.02) 


6.72 


( 26) 


2.84 


( 55) 


0.39 


07) 


0.06 (0.03) 


(0,0) 


12 C 33 S(5-4) 


242913.680 


32 


<0.06 
















(0,0) 


13 C 32 S(5-4) 


231220.996 


35 


<0.06 
















(0,0) 


HCO+(3-2) 


267577.625 


29 












15.8 


.49) 


3.56 (3.40) 


(0,0) 


H 13 CO + (3-2) 


260255.480 


29 


0.41 (.07) 


6.54 


(.08) 


1.85 


(.21) 


0.81 


.08) 


0.05 (0.05) 


(0,0) 


NH 3 (1,1) 


23694.496 


40 


3.80 (.10) 


6.17 


(.02) 


1.95 


(.05) 


25.2 


.10) 


1.30 (0.10) 


(0,0) 








2.50 (.10) 


5.78 


(.03) 


1.35 


(.09) 


10.6 


,10) 




(0,-180) 


NH 3 (2,2) 


23722.631 


40 


1.30 (.10) 


6.07 


(.07) 


2.44 


(.20) 


5.20 


.10) 




(0,0) 








0.50 (.10) 


5.95 


(.13) 


1.35 


(.16) 


0.95 


.30) 




(0,-180) 


NH 3 (3,3) 


23870.130 


40 


0.48 (.02) 


6.33 


(.07) 


3.50 


(.22) 


1.94 


.07) 




(0,0) 


NH 3 (4,4) 


24139.417 


40 


0.08 (.01) 


6.62 


(.36) 


3.77 


(.79) 


0.34 


.06) 




(0,0) 



* Columns 4 to 7 are from Gaussian fitting, except for CO(3-2), (2-1) and HCO + (3-2) for which only total integrated intensities 
are given due to self-absorption features in the spectra. Column 8 gives the 'averaged' optical depth over the line profile (see 
text). Offset positions relative to 120231 are presented in the last column. Note that column 6 gives the intrinsic full width to 
half power (FWHP) line width which has been deconvolved with the corresponding spectral resolution. Hyperfine structure has 
been taken into account fitting the C 17 0(2-l) and NH3 lines. Upper limits to the main beam brightness temperatures are 3a. 
The errors (standard deviations) of all these parameters are given in parentheses. 



4. Molecular Line Emission 

Results derived from our molecular line observations are 
summarized in Table|l], which lists: molecular line, rest 
frequency, half power beam width (HPBW) at the corre- 
sponding frequency, main beam brightness temperature, 
Local Standard of Rest (LSR) velocity, full width to half 
power (FWHP) line width, integrated line intensity, esti- 
mated optical depth averaged over the line profile, and off- 
set position relative to 1 20231. The A«i n t values presented 
in column 6 are intrinsic line widths that were deconvolved 
by the spectral resolution. The f value given in column 8 
is estimated from the ratio of integrated intensities of the 
same transition in two different isotopes, assuming isotope 
ratios of [ 13 C]:[ 12 C1=1:69.7, [ 18 0]:[ 17 0]:[ 16 0]=5. 5:1:2856 
(Wilson & Rood |1994|) , and [ 34 S]:[ 32 S]=31.3 (Chin et al. 



1996 ) for a galactocentric distance of 8.2 kpc. In the case 
of NH 3 , the peak optical depth of the main group of hy- 
perfine components is given. We assume local thermody- 
namic equilibrium (LTE) at an excitation temperature T ex 
and identical beam filling factors (unity) throughout the 
following analysis. For the line emission from linear, rigid 
rotor molecules such as CO, CS, HCO + and their isotopes, 
the beam averaged column density in cm -2 is given by 



N 



10 5 x 



3fc 2 



exp 



hvJ \ T cx + hu/6k(J+ 1) 



= / T MTi dv 



(1) 



(Scoville et al. 1986j ), where jTMsdv is the integrated 
intensity in Kkms of the J+l— >J transitions with fre- 
quency v (Hz) and optical depth f. k (ergk _1 )and h (ergs) 
denote the Boltzmann constant and the Planck constant, 
respectively, and /1 (esucm) is the permanent dipole mo- 
ment. For NH 3 column densities, see § 4.1. 



e -hu/kT B , 



4.1. NH 3 data 

Fig.|| illustrates the distribution of integrated intensity 
of the NH3 (1,1) main group of lines. The whole region 
can be divided into three subregions, with two cores lo- 
cated close to the two 870/im peaks and a ridge in be- 
tween. Such a morphology generally resembles the twin 
core system defined by Jijina et al. (1999) in their exten- 
sive study of nearby NH3 dense cores. The northern core 
has an associated IRAS source (120231), while the south- 
ern one does not. The spatial structure of the ammonia 
cores more or less follows the 870 [im dust morphology. 
With a deconvolved half-power size of 35" by 57" (0.17pc 
by 0.28 pc), the northern core is more centrally peaked 
than the southern one, which has a deconvolved size of 
about 75" (0.36 pc). To the northwest of the ridge emis- 
sion, there is a spur structure, which matches the weak 
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870 /im filament, suggesting that the continuum feature is 
real. 

While (1,1) and (2,2) line emission of ammonia is fairly 
extended in the region, (3,3) emission is more concen- 
trated around 120231 and (4,4) emission is only detected 
towards 1 20231 itself with high signal to noise ratio (S/N). 
Observed parameters of the (1,1) to (4,4) lines are listed 
in Table[| (see attached spectra in Fig.||). 

All positions detected in (1,1) and (2,2) emission 
were analysed with the usual technique (see Harju et al. 
1993| ). The rotation temperatures T ro t, the column den- 



20 



sitics iV(NH3), and the number densities ri (H2) at the 
two peak positions (0",0") and (0",-180") are calculated 
to be 15.8 and 12.9 K, 1.45xl0 15 and 7.8x 10 14 cm" 2 , 
1.9 xlO 4 and 3.5xl0 3 cm -3 , respectively. However, T rot 
derived from (1,1) and (2,2) emission usually tends to 
underestimate the gas kinetic temperature T^ n , which 
can be better approached by a T rot estimated from tran- 
sitions with higher excitation (Walmsley & Ungerechts 
1983) ; Danby et al. 1988). We thus alternatively apply a ro- 



tation diagram for a better determination of 7ki n for the 
(0",0") position where a good S/N ratio (4,4) spectrum 
was taken. Assuming a single T rot for the populations in 
all metastable levels, a best fit gives a T rot =31.8±6.2K 
from the (1,1), (2,2) and (4,4) transitions of para-NH-s. 
We can thus safely conclude a Tki n >25K for the dense 
gas within a region 40" in diameter around 120231. 

From the beam averaged H2 column density for a 40" 
area towards 120231 (see §||), we can estimate the NH 3 
fractional abundance x(NH 3 ) = [NH 3 ]/[H 2 ] -2.6xl0" 8 . 
The masses are then 43, 53 and 12 M Q for the north- 
ern core, southern core and the ridge including the spur, 
respectively. These masses agree well with those derived 
from submillimeter dust emission. The total mass of the 
whole region is therefore ~11OM0. The virial masses de- 
rived from the line width and core size are 113 and 52 M Q 
for the northern and southern core, respectively. While the 
southern core gives a virial mass very close to the masses 
derived from dust and NH3 emission, the northern core 
shows a significant inconsistency. This may be due to line 
broadening by outflowing gas which will be discussed be- 
low. 

4.2. CO data 

4.2.1. Ambient CO emission 

Spectra of selected CO isotopes taken at the position of 
120231 are shown in Fig.|. The 12 CO J = 3-2 and 2- 
1 spectra show similar line shapes, with broad-line wings 
over more than 40kms~ 1 and a deep dip at ~6 kms -1 . 
Since the off position was carefully chosen to be free of 
emission and because of the lineshapes of the 13 CO, C 18 
and C 17 profiles (see below), we believe that the dip 
feature is caused by self-absorption. The self- absorption 
between 3 and 8kms _1 , with a Tmb of about 3K and 
2K in J = 2-1 and 3-2, respectively, requires a consis- 
tent excitation temperature of T ox ~7.5K, which is likely 



10 



1 1 1 I 1 1 1 I 1 

■ 
• 

Pi 
t\ 

J \ 
f 

f \ 
J 1 

J 1 

r n 

/ nl 


1 1 l 1 1 1 l 1 

I : 
J : 


/ pi L 
y 1 


\^ 12 C0 J=3-2 " 


J j| 

/ 1 
y 

^^.^n-^^^-r"^ J 


12 C0 J=2-1 




13 C0 J=3-2 ' 






" 1 


i. \ 13 C0 J=2-1 




\ C 1B J=3-2 ' 




1 C 18 J-2-1 ' 






1 

. . . 1 ... 1 . 





-40 



-20 20 40 

Velocity (km/s) 



Fig. 3. Selected HHT spectra of CO and its isotopes to- 
wards 1 20231. The dashed vertical line marks the systemic 
velocity of WLSR=6.3kms _1 . 



the temperature of the outermost layer of the cloud. The 
13 CO, C 18 and C 17 lines are peaked at the velocity of 
the 12 CO dip (marked as a dashed vertical line in Fig. |J), 
with a Vlsr °f 6.3 km s -1 , very close to that of our NH3 
lines. From the ratio of peak temperatures of the 13 CO 
and C 18 lines, assuming a [ 13 CO]/[C 18 0] abundance ra- 



tio of 7.5 (Wilson & Rood |l99j), we find an LTE exci- 
tation temperature of 10 K for the J = 2-1 and 11.5 K 
for the J = 3-2 emission, and moderate optical depths 
of 0.96 and 0.8 at the peak velocity for the J = 2-1 and 
3-2 transitions of C 18 0, respectively. We adopt an ex- 
citation temperature of 10 K for the ambient molecular 
gas in the following calculations. Estimated from the ra- 
tio of integrated intensities of 13 CO and C 18 0, we find 
f~0.5 for the J = 3-2 and 2-1 line profiles of C 18 0. 
After applying a r-correction (multiply by a factor of 
f/(l — e _r )), we estimate LTE C 18 column densities of 
1.4xl0 16 and 7.3xl0 15 cm~ 2 for the inner 22" and 33" 
of the cloud. These lead to self-consistent fractional abun- 
dances of [C 18 O]/[H 2 ]-1.5xl0- 7 and 1.2xl0~ 7 , which are 
in good agree ment w ith the 'typical' value of 1.7xl0 -7 
(Frerking et al. |l982p . 

13 CO and C 17 were also mapped in the central part 
of the northern core (see panels a) and b) of Fig.||). C 17 
is believed to be optically thin and therefore can be used 
for mass estimation. The peak LTE column density of 
C 17 is calculated to be 8.6 xlO 14 cm 2 which gives a 
fractional C 17 abundance of 1.4xl0~ 8 . Using the in- 
tegrated emission over the mapped region, we estimate 
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Fig. 4. Contour maps of a) 13 CO(2-l), b) C 17 0(2-l), 
c) HCO+(3-2), d) H 13 CO + (3-2) and e) CS (5-4) taken 
around 120231. The contour levels are 30% to 90% by 
10% of the peak values, which are 33.3, 1.4, 17.2, 1.2 and 
6.5 K kms -1 , respectively. Line names are labeled at the 
right top corner of each panel. The plus signs mark the 
observed positions. 

a quiescent molecular cloud mass associated with C 17 
emission to be — 50M©, comparable to the mass derived 
from NH3 . Note that this mass is only a lower limit due to 
the incompleteness of the map. The deconvolved size esti- 
mated from the C 17 emission is about 40" by 75" (0.2 pc 
by 0.36 pc), slightly larger than those estimated from the 
870^m dust emission and the (1,1) line emission of NH 3 . 

4.2.2. The outflows 

Broad-line wing emission has been detected in both CO 
J = 2-1 and 3-2 transitions toward the central region in 
L870. The 12 CO wing emission, determined from a com- 
parison of C 18 and 12 CO profiles, is defined by velocities 
redder than 9 kms -1 or bluer than 3 kms -1 . Since CO 
J=2-l and 3-2 maps of the outflowing material show very 
similar bipolar structures, we present in Fig.|5| only the 
higher angular resolution CO J=3-2 data. The blue- (- 
12 to 3 kms -1 ) and the red-shifted (9 to 19 kms -1 ) lobes 
are shown in solid and dashed contours, respectively. Note 
that the integral of the rcdshiftcd wing emission was inten- 
tionally truncated to avoid contamination from a separate 
velocity component at about 20 kms -1 , which appears in 
three attached spectra at the right side of Fig.[|. 

Being centered on 120231, the outflow lobes are asym- 
metric and show some evidence of dumpiness, especially in 
the red lobe. There is some overlap between red and blue- 
shifted emission. Interestingly, close to the southern dust 
continuum peak, a pair of weak blue- and red-shifted lobe 
is also detected, indicative of ongoing star formation activ- 
ity in this region too. Velocity channel maps (Fig.^J) and 
a position velocity map along the declination axis (Fig.[?j) 
support this view. The blue-shifted lobe is, however, not 
fully covered by our map. 
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Fig. 5. Bipolar outflow of CO J = 3-2 (contours) overlaid 
on our 870/im continuum image (grey). The red-shifted 
(dashed line) emission is integrated from 9 to 19 kms -1 , 
and the blue-shifted (solid line) emission is integrated from 
-12 to 3 kms -1 . Contour levels are 30%, 40%, 90%, 
99% of the peak values, which are 35.2 (blue) and 28.8 
Kkms -1 (red), respectively. Spectra of selected positions 
(with offsets from 1 20231 in arcsec) are also shown. 

Fig.p presents a series of velocity maps in steps of 
3 kms - . We have skipped the emission between 3kms - 
and 9 kms -1 that is affected by self- absorption. As the 
maps show, the blue emission in the northern region forms 
a well-collimated lobe with 1 20231 near one end and ex- 
tends northeast, while the red emission shows two peaks, 
one at 120231 and the other 40" to the west. There are 
two clumps emerging in the 18 to 21 kms -1 channel map 
(right bottom panel in Fig. [|) , with a more compact one 
located about 2' west of 1 20231 and another one centered 
around the southern dust continuum peak SMM2. They 
are most likely cloudlets along the light-of-sight towards 
L 870, but the kinematic distances could not be estimated 
since their characteristic velocities fall into a forbidden 
zone. As shown in Fig.|7], the strong self-absorption at am- 
bient velocities in the northern core appears between two 
bright emission peaks at about 3 and 8.5 kms -1 . 

The mean optical depth in the line wings, calculated 
from 12 CO to 13 CO line ratios, is -3.0 for blue- (-12 to 
+3 kms -1 ) and —2.0 for rcdshiftcd (9 to 19 kms) wing 
emission. From the ratio of CO J = 3-2 and 2-1 bright- 
ness temperatures (the 3-2 spectra were smoothed to the 
angular resolution of the 2-1 profiles), we estimate the 
LTE excitation temperature T ox of the outflowing gas to 
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Fig. 6. Velocity channel map of CO J = 3-2. Each plot 
consists of the intensity integrated over a 3.0 km s _1 range, 
which is labeled at the right top corner of each plot. The 
lowest contour and the contour spacing are l.OKkms - . 
The positions of 1 20231 and (0" -180") are labeled as '+'. 
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Fig. 7. Position-velocity map of CO J = 3-2 emission 
along Dec. for an R.A. offset = 0. The contours and inter- 
vals are 0.3 to 2.1 by 0.6 K, and 3 to 9 by 1 K. 



by taking the characteristic flow time scale, td — Rfiow/ < 
V >, where the flow radius, i?fl ow is the lobe extension. 
Correcting the flow parameters for an inclination angle 
(to the line of sight) of 45°, we summarize the outflow pa- 
rameters in Table ||. The corresponding parameters of the 
southern outflow are estimated assuming optically thin 
wing emission, since 13 CO data are not available in this 
case. 



The parameters given in Table 2 are typical fo r low- 
mass or intermediate-mass YSOs. Cabrit & Bertout ( 1990 ) 
estimated typical errors that should approximately be 
within a factor of ~ 3 for Mfl ow , a factor of ~ 10 for F, 
and a factor of ~ 30 for £ m cch- These errors may be intro- 
duced by the uncertainties in [ 12 CO]/[H 2 ], the distance, 
line excitation temperature, inclination of the flow axis to 
the line of sight, and optical depth of the CO line. By us- 
ing intensity weighted velocities rather than the maximum 
velocities, the values of t& given in Table^j may be overes- 
timated. The flow radii Ra ow were corrected for an incli- 
nation of 45° which is likely an overestimate for the north- 
ern outflow, but an underestimate for the southern out- 
flow when we compare our observed line profiles, outflow 
structure and position- velocity maps to those of Cabrit & 
Bertout Ql99Cp . In the case of the northern outflow, effects 
caused by taking 'average' velocities and a 'standard' in- 
clination may cancel each other with respect to td- In the 
case of the southern outflow, however, they add up so that 
the td value is an upper limit for the southern outflow. We 
can therefore safely conclude that the dynamical time of 
the northern flow is longer than that of the southern flow. 



be 20 K to 40 K for an extended source and 10 K to 14 K 
for a point source. We therefore take T ex =20K for our 
outflow parameter calculations described below. 

Assuming that the CO emission is near LTE and that 
the 12 CO fractional abundance is [ 12 CO]/[H 2 ]-10" 4 , the 
outflow parameters were estimated following a method 
discussed by Cabrit & Bertout ( |1986| , |199Cp . The mass 
of the outflowing gas is estimated from the integrated CO 
J = 3-2 intensities of the red- and blueshifted lobes. The 
momentum P and the energy E are given by Mfl ow < V > 
and Mfl ow < V > 2 /2, respectively, where < V > is 
the intensity-weighted velocity. We can also get the en- 
trainment mass rate M — Mfl wAd> the force in the flow al 
F = P/td, and the mechanical luminosity L me ch — E/td, region, 



4.3. HCO+ and H id CCT 

The HCO + and H 13 CO + J = 3-2 emission lines were mea- 
sured in the central area of the northern core (see panels 
c) and d) in Fig.|J). Fig. || shows the spectra taken towards 
120231. Strong HCO + self-absorption is observed at the 
peak velocity of the H 13 CO + line (~6.5kms _1 ). Weak 
blue- and red-shifted line wing emission is seen in the J 
= 3-2 line of HCO+ that can be interpreted as additional 
evidence for outflowing gas. The asymmetry in the line 
showing enhanced blue-shifted emission characterizes all 
spectra within a radius of ^40" (~0.2pc) around 120231. 
This is indic ative of possible infall motion that will be 
discussed in S 3.4. 



The total integrated intensity map of H 13 CO + shows 
a very compact core with a deconvolved half-power size 
of about 20" by 35" (0.1 pc by 0.17pc; see panel d in 
Fig.^). With the same method used above, we estimate a 
fractional HCO + abundance of 9.6xl0~ 10 , about an order 
of magnitude lower than the typical value of ^8x 10 -9 in 



dark clouds (Irvine et al. 1987), but close to the inferred 
value of 1 x 10~ 9 for envelopes of submillimctcr continuum 
sources in the Serpens molecular cloud (Hogerheijde at 
1999). This may indicate a depletion of HCO + 



in this 

i.e. some molecules are frozen out onto the dust. 
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Table 2. Physical parameters of the two CO outflows" 



Characteristic 


-^fiow 


< V > 




Mhow 


M" 


P 


F 


E 


i'mcch 






(pc) 


(kms" 1 ) 


(ioV) 


(M Q ) 


(Moyr- 1 ) 


(Mgkms- 1 ) 


(Mokms^yr- 1 ) 


(10 45 ergs) 


(Lo) 




Blue 


0.6 


10.0 


6.0 


1.3 


2.2xl0~ b 


12.7 


2.1xl0~ 4 


1.3 


0.24 


North 


Red 


0.4 


7.8 


5.1 


0.8 


1.6xl0" 5 


6.0 


1.2xl0~ 4 


0.5 


0.10 




Total 


0.5 C 


8.9 C 


5.6 C 


2.1 


3.8xl0~ 5 


18.7 


3.3xl0" 4 


1.8 


0.34 




Blue 


0.3 d 


7.2 


4.0 


0.07 


1.8xl0 -6 


0.50 


1.3xl0~ 5 


0.04 


0.007 


South 


Red 


0.3 


7.2 


4.0 


0.08 


2.0xl0~ 6 


0.58 


1.4xl0~ 5 


0.04 


0.009 




Total 


0.3 C 


7.2 C 


4.0 C 


0.15 


3.8xl0~ 6 


1.08 


2.7xl0" 5 


0.08 


0.016 



An inclination angle of 45° has been taken into account for all parameters listed in this table. b Note the M here is the 



entrainment rate. c The average of two lobes. Due to lack of coverage in the southern outflow, Rh 
is assumed to be equal to that of its redshifted counterpart. 



for the blueshifted lobe 
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Fig. 8. J = 3-2 spectra of HCO + (dark line) and 
H 13 CO + (x5) (light line) taken at 120231. The channel 
spacing is 0.28 km s -1 . 




20 -20 

R.A. Offset (arcsec) 



The mass associated with H 13 CO is ~50M Q , consistent 
with the mass derived from C 17 0. 

As already mentioned, broad-line wings were detected 
in the HCO + J = 3-2 spectra. The outflowing emission 
shown in Fig.^ was integrated over a velocity range of 3 
to 5kms _1 for blueshifted (solid contours) emission and 
7.5 to 9kms _1 for redshifted (dashed contours) emission, 
respectively. By assuming LTE at an excitation temper- 
ature T ex =20 K and optically thin outflow emission, we 
estimate the total mass of outflowing gas traced by HCO + 
in the velocity range mentioned above to be 6 M© , almost 
three times the gas flow traced by CO. This is mostly due 
to the inclusion of the low velocity outflowing gas which 
could not directly be identified in the broader line profiles 
of CO. We thus find a total mass of 8M Q for the out- 
flow. Although such a mass estimate depends on uncertain 
molecular abundances, we can conclude that a significant 
amount of the outflowing gas is observed at low velocities 
and that the total outflowing gas mass is severely under- 
estimated by the CO data alone. 

4.4. CS and C 34 S 

In contrast to CO and HCO + our CS profiles show no 
sign of self-absorption. With a critical density of about 
7xl0 6 cm -3 , the J = 5-4 transition of CS and its iso- 



Fig. 9. Bipolar outflow (contours) overlaid on the to- 
tal integrated intensity map (grey) of HCO + J = 3-2 
taken around 120231. The blue- (solid) and the redshifted 
(dashed) emission is integrated from 3 to 5 kms -1 and 7.5 
to 9 kms -1 , respectively. Contour levels are 30%, 40%,..., 
90% of the peak values, which are 5.7 and 7.2Kkms~ 1 , 
respectively. Grey scales are 4.4 to 19.8 by steps of 
2.2Kkms _1 . 
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Fig. 10. Dual sideband spectrum with the C 34 S J = 5- 
4 line from the signal sideband and the other lines from 
the image sideband. Identified molecular transitions are 
marked. The channel spacing is 0.62 km s -1 . 
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topes traces the very dense molecular gas. While CS and 
C 34 S were detected towards 120231, C 33 S and 13 CS re- 
main undetected with an r.m.s level of ^20 mK and a 
channel separation of 0.3kms _1 . A Gaussian fit to the 
CS J = 5-4 profile at the (0",0") position gives a veloc- 
ity of 6.3 kms - , the same as that for the C 18 J = 2-1 
line but different from that of CS J = 2-1 (Bronfman et 
al. 1996). We believe that this is mainly due to different 
optical depths, i.e., optically thick emission in the J = 2- 
1 line and optically thin emission in the 5-4 line. A nine 
point CS map was also taken towards the southern core 
with the peak located at (0",-180"). As shown in Table [j], 
the CS emission lines are weaker and narrower towards the 
southern peak, and the emission appears to be much more 
compact than that in NH 3 . The integrated CS intensity 
map (Fig.^, panel e) shows a slightly elongated structure 
in N-S direction with a peak about 10" north of 120231. 

Assuming the same excitation temperature for the J 
= 5-4 lines of CS and C 34 S, we can estimate the CS op- 
tical depth for a given isotopic abundance ratio R (e.g. 
Zinch enko et al. pgi). With R = 31.3 for 120231 (Chin 



et al. 1996 ), the excitation temperature, the peak optical 
depth and the mean optical depth averaged over the line 
profile are estimated to be 6.5 K, 2.2 and 1.8. We then 
derive a CS fractional abundance of 7.8xl0 -9 , and a core 
mass of about 26 M . 




-40 -20 

Velocity (km/s) 



Fig. 11. Three 22 GHz H 2 maser spectra taken at 
1 20231 with the epoch labeled at the left top corner of each 
panel. The channel spacing is 0.13kms _1 . The full ob- 
served velocity range is ~-110 to 140 km s -1 . Only those 
parts of the spectra are shown that contain detected emis- 
sion. 



While C 34 S J = 5-4 was detected in the lower side- 
band, CH 3 OH5(l)-4(l)A- (243.9158GHz), H 2 CS7i, 6 - 
6i, 5 (244.0478GHz) and C 3 H 2 3 2 ,i-2i l2 (244.2221GHz) 
were detected in the upper sideband with peak intensi- 
ties of 0.19, 0.10 and 0.14K and integrated intensities of 
0.63, 0.57, and 0.27Kkms~\ respectively (see Fig.^). 
Since these lines are all from non-linear molecules, equa- 
tion (1) cannot be used to estimate the column densi- 
ties. Instead, we have used a relationship between the 
energy of the upper level an d the integrated line inten- 
sity (see e.g. Groesbeck et al. 1994| ) for an optically thin 
transition, adopting the line strength and partition func- 
tion information from the JPL line catalog (Pickett et al. 



1998 and online at http:/ /spec.jpl. nasa.gov ) . Under the as 



sumption of LTE and optically thin line emission, taking 
T CX =10K, we estimate CH 3 OH, H 2 CS and C 3 H 2 column 
densities to be 6.6xl0 14 , 9.0xl0 13 and l.lxlO 13 cm" 2 , 
and thus fractional abundances to be 9.7xl0~ 9 , 1.3xl0 -9 
and 1.6xl0~ 10 , respectively. These abundances are lower 
limits. If the lines are not optically thin, column densities 
and abundances would be higher. Note that abundances 
are estimated assuming that molecular lines are sampling 
the same volume as the 870 /im dust continuum emission 
(see §|3|). This does not necessarily hold, since different 
molecular species and even transitions of the same species 
have fairly different critical densities. Multi-line observa- 
tions and radiative transfer modeling are needed to test 
this hypothesis for 120231. 



4.5. Water Maser Emission 

Fig. [y] shows 22 GHz water maser spectra taken with the 
Effelsberg 100-m telescope. As shown in the spectra, the 
velocity ranges of maser emission are between -40 kins -1 
and +12kms~ , with more blue features than red ones 
relative to the systemic velocity of 6.3 km s -1 . The maser 
spectra also show time variability on a timcscalc <3 
months. While the total H 2 maser luminosity decreased 
by a factor of ~3 from 40.8Jykms~ 1 in July 1999 to 
14.3Jykms' i in July 2000, some individual maser com- 
ponents showed even stronger variations. Combining our 
data w ith the water maser spectra reported by Brand et al. 
( 1994 ), the total H 2 luminosity Lh 2 Oi assuming isotropic 
emission, varied between 3.3xlO- 7 andl.lxlO~ 6 L Q . The 
upper limit is consistent with what one expects for a 
maser associate d with an IRAS source with Lfir ,^10 2 Lq 
(Henning et al. 1992] ). From the mapped maser emission 
in April 2000, almost all the maser components, e.g. those 
at -37.1, -8.26, 4.33, 7.59 and ll.QOkms" 1 , are centered 
at the IRAS position within a pointing accuracy of ~5", 
except the 2.48kms _1 feature which peaked at (0",10"). 
We also observed the (0",-180") position (i.e. the south- 
ern NH 3 and sub- mm continuum peak), but no signal 
was seen at an rms noise level of 0.1 Jy (channnel spac- 
ing: O.lSkms^ 1 ), which yields an upper limit of Lh 2 o ~ 
3xlO- 10 L CT . 



10 



2 - 



X 



R.Q. Mao, et al.: Molecular Gas and Star Formation in Lynds870 

Table 3. 2MASS NIR properties of embedded sources 
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Fig. 12. Color-color diagram of the sources extracted from 
the 2MASS data. The circles represent YSOs. The solid 
curves show the locus of the unreddened main sequence 
and giant stars, the two parallel dashed lines form the 
reddening band and confine the region in which stars with 
normal photosphere fall. A solid circle marks the positions 
of star IRS1 (=120231). Sources not detected in J band 
are shown as open circles. Arrows indicate the lower lim- 
its, for which the J band completeness limit of 17 mag is 
adopted for J magnitudes. YSO candidates are labeled 
with sequence number from 2 and 12 according to their 
offset from 120231. In the lower right corner the average 
photometric error cross is indicated. 



5. 2MASS NIR photometry 

Stars within 8' diameter of 120231 are extracted 
from The Two Micron All Sky S urvey (2MASS, see 



http://www.ipac. caltech. edu/2mass/ for details) catalogue 



which contains J, H and K s magnitudes of stars. 
Completeness limits in magnitude are close to 17, 16 
and 15 mag in the J, H and K s bands, respectively. The 
J — H I H — K s color-color diagram presented in Fig. [l2| 
includes stars detected at least in H and K s with magni- 
tudes of 6.0< m(K s )<14.3 and photometric uncertainties 
less than 0.2 mag in all detected bands. Sources satisfy- 
ing the above criteria but not being detected in J band 
are denoted with open circles. The completeness limit of 
17 mag is adopted as a lower limit of J magnitudes for 
these sources. While the solid curve draws the locus of 
points corresponding to unreddened main sequence and gi- 
ant branch stars (Koornneef 1983), the two parallel dashed 
lines form the reddening band and are bound to the region 
of stars with normal photospheres (see Lada & Adams 



1992). The sources lying to the right side of the redden- 



ing band are believed to be mostly YSOs (Lada & Adams 



1992; Lada et al. 1993). Careful scrutiny of the region re- 



No. 


R.A.(2000) 


DEC. (2000) 


J 


H 


K, 




(h m s) 


( ° ' " ) 


(mag) 


(mag) 


(mag) 


IRS1 


20 25 07.2 


+34 50 05 


13.37 


11.21 


9.76 


2 


20 25 07.3 


+34 49 56 




15.54 


13.33 


3 


20 25 04.8 


+34 50 14 




14.71 


12.89 


4 


20 25 07.0 


+34 49 31 




15.12 


13.79 


5 


20 25 04.3 


+34 50 08 




14.49 


12.66 


6 


20 25 06.6 


+34 49 15 


_ 


15.60 


13.90 


7 


20 25 07.4 


+34 51 23 




15.58 


14.09 


8 


20 25 03.4 


+34 51 20 




15.70 


13.85 


9 


20 24 58.3 


+34 49 39 




15.76 


13.70 


10 


20 25 05.2 


+34 52 13 




15.29 


13.89 


11 


20 25 05.6 


+34 47 36 




15.54 


13.88 


12 


20 25 02.9 


+34 46 13 




15.40 


13.53 



with sequence numbers in Fig.Q see also Figj|). Since al- 
most all these YSOs except for IRS1 (solid circle) were not 
detected in J, their positions in the color-color diagram are 
still uncertain. Deep photometric observations are there- 
fore needed to confirm their YSO nature. Table || includes 
their near-infrared photometric data. From its location on 
the color-color diagram, we can estimate the extinction of 
IRS1 to be ~15, much lower than that derived from the 
peak emission of the 870/im dust continuum (Ay ~ 100). 
Although the Ay estimate using dust emission is model 
dependent (see §||), the difference in the extinction esti- 
mates is most likely beyond the uncertainties induced by 
dust properties. This may suggest that the core is struc- 
tured and/or that IRS1 is not at its center, in which case 
SMM1 and IRS1 may not refer to the same source. This 
is further supported by the fact that IRS1 is located ^8" 
(~0.04pc) northwest of the 870/xm peak (SMM1; see also 
~1. 



6. Discussion 

6.1. Physical properties of the two cores 

Dense cores are the birthplaces of stars. Their study pro- 
vides deep insights into the initial conditions of cool dense 
gas clouds that determine important stellar parameters 
like mass, metallicity, and initial rotation period. Jijina et 
al. (1999) compiled a catalog of 264 dense cores from NH3 
observations, and found that the presence of cluster asso- 
ciations has a dramatic effect on core gas and YSO prop- 
erties. The line widths, kinetic temperatures and core sizes 
for cores in clusters are much larger than the correspond- 
ing physical parameters for cores in isolated star-forming 
environments. 

The intrinsic line widths averaged over the NH3 cores 
are 1.79 and 1.20kms _1 for the northern and the south- 
ern source, respectively. These values are very close to 
the cluster associated median values of 1.20 km s -1 from 



veals 12 YSO candidates (represented by circles labeled 



the Jijina samples (samples IV and V in Jijina 199£) and 
larger than the median values for those samples with- 
out association (samples I, II and III in Jijina et al. 
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199£). Moreover, turbulence in both cores is predomi- 
nanantly non-thermal, since for the molecule of mean 
mass (=1.36hih 2 ) thermal line widths Avt (0.63 and 
0.49 km s -1 for the northern and the southern core, re- 
spectively) determined from the kinetic temperature are 
small relative to the nonthermal line widths Awnt (1-77 
and 1.18kms -1 ), obtained by subtracting in quadrature 
the thermal line width from the intrinsic line width aver- 
aged over the core. Further supported by substantial core 
masses, the two NH3 cores are reminiscent of the 'turbu- 
lent, massive dense cores' forming clusters or groups of 
stars in the Jijina et al. samples, although none of the two 
cores seems to have more than 30 embedded stars asso- 
ciated. This is supported by the distribution of the YSO 
candidates shown in Fig.[l], in which at least six (IRS1, 2 
to 6) appear to be associated with the northern core and 
another two with the southern one (11 and 12). 

Towards the northern core, there seems to be a tem- 
perature gradient from the core center to the outermost 
layer. The detection of the (4,4) line of NH3 yields a ki- 
netic temperature of ~30 K at the position of 1 20231. This 
value is an average over the beam size of 40" (~0.2 pc) and 
is likely a lower limit for the region immediately surround- 
ing 120231, since the (4,4) line was only detected towards 
the IRAS position and may not have been spatially re- 
solved. The average kinetic temperature over the rest of 
the core is ^15K. For regions where NH3 is detectable, 
the density is on the order of 10 4 cm -3 . At lower densi- 
ties (10 2-3 cm -3 ), an excitation temperature of ~ 11 K is 



traced by C 18 0, decreasing further to 
outermost layer of the clouds (see 



■ 7.5 K towards the 



4.2.1). 



6.2. The driving sources of the outflows 

The far-infrared luminosity of the IRAS source, calculated 
from the IRAS flux, is Lfir=196 L (at a distance of 

1 kpc). This luminosity is higher than the minimum value 
set by the mechanical luminosity requirement. By adopt- 
i ng th e M — Lboi relationship of Shepherd & Churchwell 
( 1996|) , the luminosity of the outflow source in 120231 
is predicted to be greater than 114L Q . Within the er- 
ror range, it suggests that the infrared source in question 
can provide sufficient mechanical power to drive the CO 
outflow. 

Source IRS1, located at (3", 0"), coincides well 
with the IRAS source 120231 with iT s =9.76mag and 
H — AT s =1.46mag. Immediately next to IRS1, a fainter 
(X s =13.33 mag) companion source 2 was found at (4", 
-11"), appearing much redder with H — K s =2.21 mag. 
Are they physically related? And which one is responsi- 
ble for the outflow? Considering the dumpiness in outflow 
structure, are perhaps both stars driving the outflow? It is 
worth noting that the position difference between source 

2 and SMM1 is 5", which is marginally within the point- 
ing uncertainty of the 870/mi observations. The question 
then arises on how IRS1, source 2 and SMM1 are physi- 
cally related. Considering the relatively large beam (22") 



of the 870/im observations, higher resolution dust con- 
tinuum and molecular line observations as well as some 
deeper NIR images are needed to answer these questions. 

In the southern core, we did not find such YSO can- 
didates close to SMM2. Although there is a much weaker 
(14.13 mag) K s source located at (-3", -180"), its nature 
is unknow since it was not detected in J and H and we 
cannot place it on the color-color diagram. Most likely 
the driving source of the southern outflow is still deeply 
embedded in the cloud. It is also less luminous than its 
counterpart (s) of the northern outflow. 

6.3. Evolutionary status of the two cores 

Protostars evolve from Class to Class I, II, III. Class 
sources, the youngest protostars, are characterized by 
M env /Lboi > O.1M /L Q (M env is the envelope mass), no 
detectable emission at wavelengths shorter than ~ 10 /im, 
and the presence of outflow (Andre et al. 



1993, 2000). At 



present, 42 Class sources have been identified (Andre et 
al. p00| ). 

A reasonably good estimate of M onv can be provided 
by the strength of the optically thin 870 /im continuum 
emission in the 22" beam (~0.1pc, the typical size of 
protostellar envelopes in the Taurus cloud) . Following the 
method used in §pl we can estimate M env for SMM1 and 
SMM2 to be about 9 and 14 M , with their peak flux den- 
sities of 1.67 and 0.60 Jy beam -1 , respectively. It is, how- 
ever, worth noting that (3 = 1.5 (Ks7o M m = 0.02 cm 2 g -1 
correspondingly) , a value lower than the one used for mass 
calculations in §^|, is found to be more suitable to esti- 
mate the m asses traced by the dust around protostars (e.g. 
Visser et al. |l99?j , 2001 ). Moreover, the averaged dust tem- 
perature within a 22" beam could be higher than the ones 
used in § H Both effects lead to lower masses. Considering 
additionally that 120231 may not be at the core center, 
we can safely conclude that M onv /Lboi < 0.04Mq/L q for 
SMM1, typical of a Class I source. The detection of 1 20231 
at 12 /im and 2 /im and the moderate strength of its CO 
outflow also suggest SMM1 to be a ~ 200 L Class I YSO. 

To better constrain the properties of SMM2, we first 
estimate the lower limit to M env . «870/im = 0.03 cm 2 g -1 is 
most likely an upper limit if we take the assumption from 
Hildebrand ( |1983| ), i.e. k x = 0.1(250/A(/im)) /3 cm 2 g -1 , 
with f3 = 1, which is the lower limit in star forming re- 
gions at the millimeter and submillimcter wavelengths. A 
dust temperature of 30 K is most likely an upper limit 
to SMM2. These parameters yield a lower limit to the 
envelope mass of 1M . Meanwhile, since there is no de- 
tectable IRAS source associated with SMM2, the sensitiv- 
ity of the IRAS instruments (0.5 Jy at 12, 25, and 60 /im, 
1.5 Jy at 100 /im) allow us to estimate an upper limit of 
^boi ~6L/0 for the luminosity of the possible source(s) 
embedded in the southern core. The M env /Lboi ratio for 
SMM2 is therefore comfortably greater than O.IMq/Lq, 
typical of a Class source. Combined with the fact that 
SMM2 also harbors a bipolar outflow, the southern core 
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is more likely to contain a Class object with 1 M Q 
M onv ^14M Q . 

Since outflows are found in both cores, we can also 
compare the evolutionary phases of the two cores by plac- 
ing these two sources on the outflow efficiency versus 
Menv/^bo! diagram (see fig. 7 of Bontemps et al. 1996| ). 



The outflow efficiency is the outflow force divided by the 
radiative momentum flux (Lboi/c) and is ^80 (or ^380 
also accounting for the low velocity flow traced by HCO + ) 
for the northern outflow and >220 for the southern one. 
Accounting for the different M env /L^ values of the two 
cores, the northern one lies in the transit region between 
Class and I, while the southern one lies distinctly in 
the region occupied by Class sources. This, in agreement 
with the fact that the dynamical time of the northern flow 
is longer than that of the southern flow (see § 4.2.2), fur- 
ther supports the protostellar status of SMM1 (120231) 
and SMM2. 

We thus conclude that the southern core is in an earlier 
evolutionary stage than the northern core. One may ask 
why there is no water maser detection in the southern core, 
since Class sources are believed to be favorable sites for 



such masers (Furuya et al. 2001). An answer could be that 
the maser is too weak to be detected at an rms level of 
0.1 Jy, since according to the £h 2 — -^boi correlation found 



by Furuya et al. (2001), we expect for a protostar with 
-^boi < 6Lq a luminosity of Lh 2 o Si 1.7 x 10~ n Lq. This is 
about two orde rs of magnitude below the estimated upper 
limit (see § |4.5| ). On the other hand, the non-detection of 
a H2O maser in the southern core is not inconsistent with 



the recent statistical study by Zhang et al. (|200l[ ) who 
found that the outflow may develop before the appearance 
of a water maser although their study is based on a sample 
of massive star forming regions. It is therefore also possible 
that the southern core has not yet reached an evolutionary 
stage to form a luminoous H2O maser, while an outflow 
was already activated. 

The above disscussion is based on the assumption that 
SMM1 and IRSI refer to the same source as 120231, but 
this is not necessarily the case. As we have mentioned in 
§H and § |6.2j , the differences in extinction (A v ~ 100 versus 
~15) and in position between SMM1 and IRSI (-8" or 
0.04 pc) raise the possibility that the two sources are dis- 
tinct objects of the northern "cluster" . In such a case, Lboi 
of SMM1 would be significantly lower than the value esti- 
mated from FIR flux densities of 120231, and as a result 
SMM1 might still be a Class object and driving source 
of the northern outflow. High angular resolution observa- 
tions would be required to draw definite conclusions. 

6.4. Velocity gradients and the asymmetric line profiles 

There exists a velocity gradient from north to south in all 
molecular spectra, indicating an internal motion. In the 
case of the NH3 lines, the velocity increases monotonously 
from 5.2kms _1 at (0",80") to 6.8kms _1 at (0",-80"), 
where there seems to be a turning point. Further south 





Fig. 13. Position-velocity maps of 13 CO J = 2-1 (left) 
and NH3 (1,1) (right) emission along Dec. for R.A. offset 
= 0. The contours are 10% to 90% by 10% of the peak 
emission (9.2 K and 4.1 K for 13 CO J = 2-1 and NH 3 (1,1), 
respectively). 

the velocity shifts back to and stays almost constant at 
6.0kms _1 . Such a gradient is also found in 13 CO J = 2- 
1, where the velocity shifts from 5kms _1 at (0",75") to 
6.3kms _1 at (0",-45"). The CS J = 5-4 line profile at 
the central position is quite similar to that of the 13 CO 
line, with a peak slightly skewed to low velocity. Again, 
the velocity increase from 4.5 kms" 1 at (0",75") to 6.5 
kms -1 at (0",-30"), and drops to 6.0 kms 1 further south 
just like ammonia. This indicates that a boundary ex- 
ists near (0",-100"), separating the northern cloud (with a 
'systemic' velocity ~ 6.3 kms 1 ) from the southern cloud 
(with ~ 6.0 km s^ 1 ). The gradients across the northern 
core are 2 to 4 kms -1 pc -1 , slightly larger than the typical 
rotation induced gradients of lkms -1 pc -1 for molecular 
cloud cores (Goodman et al. 1992), but within the range 



found by Jijina et al. (1999) in their sample. As examples, 
we present in Fig.[l3| the position- velocity maps of 13 CO 
J = 2-1 (left) and NH 3 (1,1) (right). Morphologically, the 
diagrams are consistent with a large cloud that shows over- 
all large scale rotation. However, such an effect could also 
be caused by two independent cores with slightly different 
velocities. 

The HCO + J = 3-2 profiles are of particular inter- 
est, owing to the presence of strong blue asymmetry (sec 
§ 13), indic ating possible infall m otion (Leung & Brown 
1977t Zhou |l992j ; Gregersen et al. |l997[ ). The spatial ex- 
tent (~ 0.2 pc) of such an asy mmetr y is not inconsistent 
with the finding of Myers et al. 200C , who detected infall 
motion over a linear scale of ~ 0.1 pc in a large number of 
starless cores and Class 0/1 protostars. However, it is hard 
to study the detailed infall structure with current spa- 
tial resolution. Contamination by rotation and outflows 
make a clearcut interpretation even more difficult. Further 
high spatial resolution observations are imperative to ver- 
ify such a possible infall. 

7. Conclusions 

We have presented detailed molecular line and sub-mm 
dust observations of the molecular cloud Lynds 870, a pre- 
viously poorly studied region. The main results are: 
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1. Substantial molecular emission has been detected 
from Lynds870. Two main cores, separated by 3', are 
found both in dust (SMM1 and SMM2) and molecular 
line emission with the northern core (SMM1) being partic- 
ularly rich in molecular lines. The continuum and spectro- 
scopic observations indicate that the northern core, cen- 
tered at IRAS 20231+3440, is warmer and denser than the 
southern core (SMM2). A total mass of ~ 70-110 M is 
estimated from the molecular line and the submillimeter 
dust emission. 

2. Outflow activities are detected in both the north- 
ern and the southern core, indicative of on-going star for- 
mation in both regions. While IRAS 20231+3440 (IRS1 
in K s ) is found to be a Class I YSO candidate and most 
likely the exciting star of the northern outflow, a more 
deeply embedded companion (source 2 in K s ) found next 
to IRAS 20231+3440 may also be a possible driving en- 
gine. The driving source of the southern outflow, still 
deeply embedded in the cloud, is most likely a Class 
YSO candidate. The southern core (SMM2) is therefore 
still in an earlier evolutionary stage than the northern 
core (SMM1). 

3. A significant amount of low velocity outflowing gas is 
hidden in the relatively broad systemic CO profiles tracing 
the quiescent gas. The low velocity outflow is, however, 
revealed by the HCO + emission. The total outflowing gas 
mass of the northern core is derived to be 8 M Q . 

4. Possible infall motions in the northern dense core 
are indicated by self-absorbed HCO + J = 3-2 line pro- 
files with a strong blue asymmetry; its spatial extent is 
-0.2pc. 

5. In a diameter of 8' around IRAS 20231+3440, there 
are in total 12 K s sources with a strong IR excess, typical 
of YSOs still embedded in their parental dust/molecular 
cloud. Combined with the outflows and H 2 masers found 
in the region, we conclude that the cloud is presently ac- 
tively forming a group or a small cluster of stars. 
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